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FOREWORD 


This  report,  which  presents  the  experimental  and  theoretical  results  of  a 
program  conducted  to  Investigate  the  supersonic  flutter  characteristics  of  a 

l 

swept-back  all-movable  surface,  was  prepared  by  the  Aeroelastic  and  Structures 
Research  Laboratory,  Massachusetts  Institute  of  Teclmology,  Cambridge  39, 
Massachusetts  for  the  Aircraft  Laboratory,  Wright  Air  Development  Center, 
Wright-Patterson  Air  Force  Base,  Ohio.  The  work  was  performed  at  MIT 
under  the  direction  of  Professor  R.  L.  Halfman,  and  the  project  was  supervised 
by  Mr.  G.  W.  Asher.  The  research  and  development  work  was  accomplished  under 
Air  Force  Contract  No.  AF  33(616) -2751,  Project  No.  1370  (UNCLASSIFIED 
TITLE)  Aeroelastic  tty,  Vibration  and  Noise,”  and  Task  No.  13479, 
(UNCLASSIFIED  TITLE)  "Investigation  of  Flutter  Characteristics  of  All-Movable 
Tails, "  with  Mr.  Niles  R.  Hoffman  of  the  Dynamics  Branch,  Aircraft  Laboratory, 
WADC  as  task  engineer.  This  research  was  started  in  January  1955  and  completed 
in  September  1956.  Additional  supersonic  flutter  testing  of  swept  all-movable 
stabilizers  may  be  performed  at  a  later  date  to  obtain  further  information. 

The  authors  are  indebted  to  Mr.  O.  Wallin  and  Mr.  C.  Fall  for  their  help  in 
building  the  models  and  in  keeping  the  experimental  equipment  in  good  order,  and 
to  Mr.  G.  M.  Falla  for  his  help  in  making  the  high  speed  photographs.  The 
authors  are  also  Indebted  to  Messrs.  A.  Heller,  Jr. ,  J.  R.  Friery  and  H.  Moser 
for  their  help  in  preparing  the  necessary  calculations,  tables,  and  figures  for 
this  report. 

Portions  of  this  document  are  classified  CONFIDENTIAL  since  the  data 
revealed  can  be  employed  to  establish  design  criteria  for  the  prevention  of  flutter 
of  swept-back  all-movable  taiis  of  aircraft  in  the  supersonic  speed  range. 
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ABSTRACT 


^  This  report  describes  the  flutter  testing  at  supersonic  speeds  of  a  series  of 
swept-back  all-moving  stabilizers.  An  attempt  was  made  to  define  the  flatter 
boundaries,  for  one  location  of  the  pitching  axis,  over  the  Mach  number  range  of 
1. 3  to  2. 1,  by  testing  at  a  number  of  different  levels  of  stabilizer  stiffness,  and 
at  a  number  of  different  pitching  frequencies. 

The  results  indicate  that  large  Increases  In  the  region  of  instability  can 
occur  due  to  the  Introduction  of  the  pitching  degrees  of  freedom.  The  test  results 
follow  the  trends  of  theoretical  calculations,  but  the  quantitative  correlation  be¬ 
tween  the  theoretical  and  the  experimental  results  is  only  fair. 
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SECTION  I 
INTRODUCTION 

This  report  covers  the  experimental  flutter  tests  and  associated  theoretical 
calculations  made  on  a  swept,  all-moving  horizontal  stabilizer  at  supersonic 
speeds.  The  configuration  tested  is  becoming  a  comrr><;n  one  for  high  speed  air¬ 
craft  and  missiles. 

At  present,  the  methods  of  theoretical  supersonic  flutter  analysis  using 
two-dimensional  aerodynamic  forces  derived  from  linearized  theory  do  not 
appear  adequate  to  predict  the  absolute  levels  of  the  flutter  boundaries. 

Reference  3  shows  that  even  for  the  simple  cantilever  straight  wings  analyzed 
in  that  report  such  analyses  give  results  that  are  conservative  in  one  Mach 
number  range  and  unconservative  in  another.  It  may  be  suspected  that  the  poor 
correlation  between  theoretical  and  test  results  shown  in  Ref.  3  arises  from 
the  use  of  two-dimensional  aerodynamic  forces  on  a  three-dimensional  lifting 
surface,  and  so  the  use  of  more  powerful  methods  of  analysis,  such  as  the 
aerodynamic  influence  coefficient  methods  of  Refs.  5  and  6,  may  improve  the 
correlation.  However,  correlations  between  theoretical  and  experimental 
results,  where  the  theoretical  calculations  have  been  based  on  three-dimensional 
aerodynamic  forces,  are  not  common  in  the  supersonic  regime.  Until  such 
correlations  have  been  made,  it  is  not  certain  that  the  added  labor  of  the  in¬ 
fluence  coefficient  methods  will  be  worth  while  in  terms  of  improved  results. 

The  designer  will  probably  rely  on  the  simpler  two-dimensional  calculation 
to  supply  a  description  of  the  trends  to  be  expected  when  various  parameters 
are  changed,  and  will  probably  depend  for  some  time  on  what  experimental  data 
is  available  or  can  be  obtained  to  define  the  absolute  levels  of  the  flutter 
boundaries. 

The  present  program  is  intended  to  define  experimentally  the  level  of  the 
flutter  boundaries  for  an  all -moving,  swept  horizontal  stabilizer.  The  canti¬ 
lever,  or  "pitch  locked,  "  boundary  is  defined  by  tests  of  the  cantilever  con¬ 
figuration  of  the  model  shown  in  Fig.  14a,  and  through  the  use  of  data  from 
previous  flutter  tests.  Various  levels  of  wing  and  pitching  restraint  stiffness 
are  then  combined  in  an  attempt  to  define  the  effect  of  the  pitching  degree  of 

Manuscript  released  by  the  authors  September  1957  for  publication  as  a 
WADC  Technical  Report. 
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freedom  on  the  flutter  boundaries,  over  the  Mach  number  range  1.  3  to  2.  0.  The 
results  are  discussed  In  Section  Q  of  this  report,  and  a  complete  compilation 
of  the  experimental  data  is  found  in  Appendix  II. 

Along  with  the  experimental  program,  a  large  number  of  theoretical 
calculations  have  been  made  on  the  basis  of  two-dimensional  aerodynamic 
coefficients,  both  supersonic  and  incompressible.  The  major  effort  was  ex¬ 
pended  on  three-degree-of-freedom  calculations  employing  assumed  wing  bend¬ 
ing  and  wing  torsion  structural  modes  and  a  rigid  pitching  mode.  Four-degree- 
of-freedom  calculations  were  also  made  which  included  an  assumed  second  bend¬ 
ing  mode  as  well  as  the  previously  mentioned  modes.  The  results  of  the  calcula¬ 
tions  are  discussed  in  Section  II,  and  the  equations  used  for  setting  up  the  cal¬ 
culations  are  described  in  Appendix  I. 
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SECTION  n 


DISCUSSION  OF  RESULTS 


I  Discussion  of  Theoretical  Results 


A  considerable  number  of  calculations  were  made  during  the  course  of  the 
program  for  the  swept  stabilizer  models  to  determine,  if  possible,  what  trends 
might  be  expected  In  the  flutter  boundaries  for  models  of  various  stiffness  levels 
and  varying  pitch  frequencies.  The  calculations  used  the  velocity  component 
method  of  Ref.  8  with  two  dimensional  aerodynamic  coefficients.  The  model  for 
calculation  of  the  aerodynamic  integrals  was  assumed  to  be  untapered  In  order 
to  avoid  variations  of  reduced  frequency,  k  =  buij/v,  along  the  span,  but  the 
mass  and  stiffness  parameters  were  assumed  to  vary  In  the  same  manner  as  the 
experimental  model.  A  description  of  the  calculations  is  found  In  Appendix  I. 


Most  of  the  theoretical  effort  was  expended  on  three-degree-of-freedom 
calculations  employing  supersonic  aerodynamic  coefficients.  The  three  degrees 
of  freedom  used  for  this  analysis  were  wing  first  bending  (parabolic).wing  first 
torsion  (linear), and  rigid  pitch  about  the  rotation  axis.  The  results  of  these 
calculations  are  given  In  Fig.  1  for  various  values  of  the  pitching  to  torsion 


2 

frequency  ratio  (  )  and  pitching  to  first  bending  frequency  ratio 

,  /  ,2  . ^ 


(<*>^  /w 
and  (t> 


In  all  of  the  calculations 


( <*>.  /w 

hj  Oj 


)  =  0.  25,  where  w.  ,  w 


1  “I  UI  '*1  -l 

refer  to  the  frequencies  in  the  assumed  uncoupled  modes.  Note  that 


the  reference  semichord  for  the  calculations,  br,  is  that  of  the  75  %  span 
station,  bp  The  reference  axes  for  measuring  the  semichord  as  well 


as  other  similar  quantities  are  aligned  with  the  stream  unless  a  subscripts 
is  used.  In  that  case  the  reference  axes  are  swept  with  the  elastic  axis  (See 
Fig.  8>.  Many  of  the  parameters,  such  as  p,  are  constant  along  the  span. 


Perhaps  the  most  interesting  feature  shown  in  Fig.  1  is  the  sharp  increase 
in  the  region  of  instability  that  occurs  below  Mach  numbers  of  about  1 . 8  for  a 
value  of  )  =  0.  20.  For  this  value  of  (w^/ui^  )^,  the  boundary  actually 

crosses  the  M  -  1.7  line  three  times  giving  the  shape  shown.  As  noted  In  Ap¬ 
pendix  I,  the  four -degree -of -freedom  calculations  Indicated  that  this  sharp  drop 
or  •bucket"  in  the  boundary  probably  occurs  because  of  a  change  in  the  mode  of 
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F ig.  1.  Flutter  porometers  V(  ■»,  Oj  77  ond  Vf/wfbp  versus  Moch  number  from  tfiree-degree-of-freedom 

supersonic  calculations. 
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flutter.  Figure  13  of  Appendix  I  shows  that  the  "bucket"  may  be  very  much  af¬ 
fected  by  the  level  of  structural  damping,  g.  to  the  extent  that  for  g  about  0.6  the 
boundary  for  (u>^/u>a  )2  =  0.20  may  follow  smoothly  the  trend  established  by  the 

1  2 

cantilever,  or  "locked"  case  (u>,  /u>  )  =  *. 

©  Oj 


For  high  pitch  frequencies  the  velocity  -  Mach  number  trend  is  similar  to 
the  locked  pitch  case  with  only  a  moderate  lowering  of  the  flutter  velocity.  For 
low  pitch  frequencies  the  velocity  -  Mach  number  trend  appears  to  be  nearly  a 
constant  equivalent  airspeed  (V’o-’V  =  constant)  through  a  wide  range  of 
Mach  number. 


For  M  1.5,  corresponding  to  a  cross  flow  Mach  number  perpendicular  to 
the  40% chord  line,  M^j,  of  10/9,  a  sharp  Increase  in  the  region  of  instability 
occurs  even  for  the  cantilever  case,  (w  /  )2  =r  ®.  This  increase  in  the 

region  of  instability  may  arise  from  the  use  of  the  linearized  supersonic  aero¬ 
dynamic  theory  at  such  a  low  cross  flow  Mach  number.  It  must  be  remembered 
that  the  Mach  number  used  for  the  determination  of  the  aerodynamic  coefficients 
of  the  calculation  is  the  cross  flow  Mach  number,  not  the  free  stream  Mach 
number. 

The  cantilever  curves,  (u^/wQ^)2  =ec,  of  Fig.  1  were  determined  from  twu- 

degree-of-f reedom  calculations  in  which  first  bending  and  first  torsion  modes 
were  used  without  the  pitching  degree  of  freedom. 

Figure  2  shows  curves  of  the  flutter  parameters  V,/u>  b-  and 

I  Ct j  U.  75 

vf/u)f  b0  75  versus  the  frequency  ratios  /xQ  1  and  (wf/u>h  )  calculated  by- 

using  incompressible  aerodynamic  coefficients  and  three  degrees  of  freedom: 
wing  first  bending,  wing  first  torsion,  and  rigid  pitch.  A  sharp  increase  in  the 
region  of  instability  for  values  of  )  less  than  about  0.  5  can  be  seen. 

More  significantly,  the  decrease  In  stability  appears  to  be  related  to  the 
near  equality  of  pitch  and  bending  frequencies.  This  effect  has  been  observed 
by  other  investigators  and  depends  on  pitch  axis  location. 

Calculations  were  also  made  using  an  assumed  second  bending  mode  along 
with  the  first  bending,  first  torsion,  and  pitch  modes  and  the  results  are  dis¬ 
cussed  in  Appendix  I.  The  addition  of  the  second  bending  mode  does  not  affect 
the  shape  of  the  flutter  boundaries  significantly  in  the  Mach  number  range  studied. 
Changing  the  ratio  of  second  bending  to  first  torsion  frequency,  (u>h  ), 

from  slightly  greater  than  1.0  to  slightly  less  than  1.0  also  has  little  effect  on 
the  flutter  boundaries.  It  appears,  then,  that  sufficient  accuracy  was  obtained 
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in  the  calculations  with  three  degrees  of  freedom  for  the  wings  studied.  As  noted 
in  Appendix  I,  the  V-g  solutions  for  the  four-degree-of-freedom  calculations  did 
furnish  valuable  insights  into  the  modes  of  flutter  and  the  effect  of  structural 

damping. 

2  Discussion  of  Experimental  Results  . 

During  the  test  program  nine  cases  of  flutter  occurred  for  the  sixteen 
configurations  tested.  Two  models  fluttered  in  a  cantilever,  or  "pitch  locked." 
condition  and  the  remaining  models  at  various  levels  of  wing  stiffness  and 
pitching  frequency. 

Reference  1  describes  the  M.  I.  T.  -WADC  supersonic  variable  Mach  number 
Blow-Down  Wind  Tunnel  facility  in  which  the  tests  were  conducted.  Reference  2 
describes  the  techniques  of  testing  that  were  used  to  obtain  the  data.  No  major 
changes  were  necessary  in  either  the  wind  tunnel  facility  or  the  testing  techniques 
to  obtain  the  experimental  data  presented  in  this  report. 

The  planform  of  the  stabilizer  models  tested  is  shown  in  Fig.  14  of 
Appendix  II.  They  incorporated  a  pitching  degree  of  freedom  with  a  pitch  axis 
perpendicular  to  the  root  chord,  64.  3%  of  the  root  chord  aft  of  the  leading  edge. 
The  stiffness  of  the  pitching  restraint  could  be  varied  at  will.  The  model  con¬ 
struction  was  similar  to  that  described  in  Ref.  2  with  a  single  spar  providing  the 
required  stiffness.  Balsa  fairings  glued  to  the  spar  gave  the  required  6%  thick 
double  wedge  airfoil  shape  and  suitably  spaced  lead  weights  provided  the  required 
mass  parameters.  A  more  complete  description  of  the  models  is  given  in 
Appendix  n. 

Before  flutter  testing,  each  model  was  given  vibration  and  static  tests.  The 
results  of  these  tests,  as  well  as  the  tabulated  results  of  the  flutter  tests,  are 
contained  in  Appendix  n.  With  the  pitching  mechanism  "locked  out,  "  the 
cantilever  condition,  the  lowest  natural  modes  of  vibration  were  determined  lor 
each  model.  In  general  three  modes  were  easily  excited,  the  first  bending, 
first  torsion,  and  second  bending  modes.  The  first  bending  mode  and  first 
torsion  mode  determined  in  this  maimer  were  used  to  plot  the  flutter  data  of 

Figs.  3  and  4  are  the  and  of  the  figures.  The  rigidities  in  bending  and 
N  N 

torsion,  EIr  and  GJr,  at  the  root  were  also  determined  for  most  of  the  models  in 
the  cantilever  condition.  This  data  is  not  too  satisfactory  since  it  is  difficult  to 
assess  accurately  the  effects  of  root  fitting  deformation.  As  can  be  seen  from 
Table  2  there  seems  to  be  considerable  scatter  in  the  EIr  and  GJ  data  since 
models  with  essentially  the  same  cantilever  frequencies  appear  to  have  widely 
different  values  of  EIr  and  GJr.  With  the  pitching  mechanism  in  operation 
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vibration  data  was  also  taken  for  various  pitch  restraint  stiffnesses.  In  general, 
only  the  first  three  modes  of  vibration  could  be  excited  easily  as  can  be  seen 
from  the  data  of  Table  5.  This  data  furnished  the  coupled  vibration  frequencies 
and  for  the  plots  In  Figs.  6  and  7.  Influence  coefficient  data  was  also 
taken  with  the  pitching  mechanism  in  operation.  The  uncoupled  pitch  frequency, 
Wq,  was  determined  from  the  measured  rigidity  of  the  pitch  mechanism  and 
from  the  measured  total  mass  moment  of  inertia  of  the  wing  and  root  fitting,  1^ . 
A  few  of  the  frequencies  so  determined  were  checked  by  fitting  a  rigid  disc  of 
known  moment  of  inertia  to  the  flexure  and  measuring  the  resulting  vibration  fre¬ 
quency.  The  check  on  frequencies  was  satisfactory.  Pitching  frequency  data  can 
be  found  in  Table  3  of  Appendix  II,  while  the  frequency  data  and  all  of  the  flutter 
data  is  summarized  in  Table  4. 


Figures  3.  4,  and  5  compare  the  experimei  ■]  flutter  data  and  the  theoretical 

predictions  when  plotted  versus  Mach  number,  .t  is  presumed  that  w  ,  the 

°N 

first  measured  cantilever  torsion  frequency,  corresponds  fairly  closely  to  the 
uncoupled  first  torsion  frequency,  u>  ,  used  as  a  parameter  in  the  calculations 


and  similarly  that  w.  corresponds  closely  to  w.  .  Since  the  different  models 
nN  hl 

fluttered  at  somewhat  different  relative  densities,  n{,  and  since  the  value  of  p 
used  in  the  theory  is  lower  than  for  most  experimental  points,  the  factor 
l/Vpj  has  been  included  in  the  ordinates  to  reduce  the  effects  of  these  variations. 


The  tests  of  the  SWS-2  model,  which  fluttered  in  a  locked  configuration, 

along  with  the  data  of  Ref.  3  were  used  to  establish  the  cantilever,  or  bending- 

torsion  flutter  boundary;  (  -u  /uv  )2  or  (w  /u»  )2  =  ®.  (The  SWS-ld  model 

v  "N  v  “n 

also  fluttered  in  a  cantilever  condition  but,  since  the  vibration  data  of  Table  5 
shown  that  this  model  had  a  low  torsion  frequency  quite  different  from  the 
rest  of  the  stabilizer  models,  it  was  used  only  as  a  guide  in  drawing  the  "locked” 
boundary. ) 


The  SWS-1  series  of  models,  had  a  slightly  higher  stiffness  level  than  SWS-2 
and  thus  had  a  margin  of  safety  of  about  7%  in  bending -torsion  flutter.  The 
margin  of  safety  Is  defined  as  the  ratio  of  &  necessary  to  prevent  flutter  in 

the  cantilever  condition  to  the  of  the  actual  model.  The  SWS-1  series  models 

N 

were  flutter-free  in  the  cantilever  condition  but  when  the  pitch  frequency  was  low¬ 
ered  to  about  98  cps,  flutter  occurred  at  M  --  1.35,  as  can  be  seen  from  the 
SWS-1 -98  model  test  point.  Two  other  SWS-1  series  models  were  flutter  tested 
at  lower  values  of  pitch  frequency,  the  SWS-lc-48  and  the  SWS-le-74  models. 

The  vibration  data  shows  that  these  models  were  similar  to  the  SWS-1 -98 


vVADC  TR  50-285 


8 


CONFIDENTIAL 


CONFIDENTIAL 


CONFIDENTIAL 


Flutter  pedometer  V^/ ^0.75^/  0f'v*r*us  Moch  number  from  experimental  tests  ond 

with  theory. 
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model  In  a  cantilever  condition.  The  SWS-lc-48  model  fluttered  on  injection 
when  practically  in  the  tunnel  whereas  the  SWS-le-74  model  fluttered  in  the 
middle  of  a  test  run. 


The  three  flutter  points  for  the  SWS-1  series  models  cover  quite  well  the 
Mach  number  range  available  in  the  wind  tunnel  so  that  further  tests  of  this 
series  of  models  at  intermediate  values  of  the  pitching  frequency  were  not 
attempted.  Instead  a  third  series  of  models,  SWS-3,  were  designed  with  a 
margin  of  safety  In  bending-torsion  flutter  of  about  45% based  on  the  curves 
of  Figs.  3,  4  and  5.  In  order  to  achieve  the  higher  frequency  and  stiffness 
level  required  by  the  increased  margin  of  safety  without  increasing  the  model 
thickness  ratio,  it  was  necessary  to  modify  the  design  parameters  of  the  SWS-3 

models.  The  frequency  ratio  (wK  /  ui  )  was  lowered  from  an  average  of  0.  29 

nN  °N 

to  0.  26  rather  than  change  the  mass  parameters.  The  test  data  and  calculations 
of  Ref.  3  show  that  there  is  little  variation  in  the  level  of  the  cantilever  flutter 

boundaries  for  straight  and  swept  wings  for  variations  in  (w.  / w  )  over  this 

rN  “N 

range. 


Three  of  the  SWS-3  series  models,  SWS-3b-53,  SWS-3a-63,  and  SWS-3c-74 
fluttered  on  or  very  close  to  injection.  The  SWS-3c-74  fluttered  when  fully  in 
the  tunnel  but  before  the  Mach  number  had  started  to  change  and,  therefore,  is 
not  shown  as  an  Injection  flutter.  The  SWS-3a-63  and  the  SWS-3b-53  were  almost 
In  the  tunnel  when  flutter  occurred  and  are  shown  as  injection  flutter.  In  sketch¬ 
ing  the  experimental  boundaries, the  data  for  SWS-3d-87  and  SWS-3c-74  were 
relied  on  more  heavily  than  the  data  for  SWS-3a-63  and  SWS-3b-53. 


The  SWS-3  series  data  show  that  there  can  be  a  very  large  increase 
in  the  region  of  Instability  if  the  ratio  is  near  unity.  In  fact,  it 

2  ^ 

appears  that  for  a  given  value  of  (ui./w.,  )  the  stiffer  SWS-3  models  will  flutter 

0  «N 


at  higher  Mach  number  than  their  SWS-1  counter  parts.  Thus,  the  experimental 

2 

boundaries  for  a  given  ( u>.  /u>.  )  appear  to  bend  back  and  form  deep  "buckets" 

o  nN 

in  the  curves  just  as  they  do  for  the  calculated  results.  Ingeneral,  however, 
the  calculations  predict  larger  regions  of  instability  than  the  experimental 
results  indicate. 


It  is  interesting  to  note  that,  the  SWS-3  series  flutter  apparently  occur  in 
a  different  flutter  mode  than  the  SWS-1  series.  Figures  18  and  19 
show  the  analysis  of  the  high  speed  movies  for  the  SWS-1 -98  and  the  SWS-3d-87 
models,  taken  from  the  excerpts  from  the  high  speed  movies  shown  in  Figs.  16 
and  17.  The  SWS-1 -98  model  should  have  a  different  mode  of  flutter  than  the 
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SWS-3d-87  if  the  results  of  Appendix  I  are  correct  in  that  the  "buckets”  of  the 
boundaries  are  formed  by  a  new  flutter  mode.  Examination  of  Figs.  18  and  19 
shows  that  while  the  relation  between  the  tip  vertical  translation  amplitude  to 
pitch  amplitude  is  of  the  same  order  of  magnitude  for  the  two  models,  the 
relationship  between  the  tip  angle  of  attack,  amplitude  and  the  pitch  amplitude  is 
much  different.  The  SWS-1-98  model  shows  a  much  larger  ratio  of  tip  angle  of 
attack  amplitude  to  pitch  amplitude  than  does  the  SWS-3d-87.  This  fact  Indicates 
that  the  flutter  mode  for  the  SWS-1-98  is  composed  of  Important  pitch -bending  - 
torsion  motions  while  the  flutter  mode  for  the  SWS-3d-87  is  mainly  pitch-bending. 

It  would  appear,  then,  that  for  small  margins  of  stability  in  bending-torsion 
flutter  the  addition  of  a  high  frequency  pitch  degree  of  freedom  causes  a  decrease 
in  what  is  essentially  a  bending-torsion  flutter  speed  largely  because  of  the 
decrease  In  the  coupled  torsion  'requency.  However,  if  (w0/a>h)  is  low  enough 
to  be  near  unity  a  bending  pitch  mode  develops  which  may  increase  the  region 
of  instability  to  as  high  as  M  -  2. 

Before  discussing  some  of  the  other  curves  drawn  from  the  test  data,  some 
attention  should  be  given  to  the  SWS-3-53  model.  This  model,  although  practi¬ 
cally  identical  with  the  SWS-3b-53  model  insofar  as  vibration  frequencies  are 
concerned,  was  tested  in  the  same  range  of  Mach  number  and  density  as  the 
SWS-3b-53  model  but  failed  to  flutter.  However,  the  structural  damping  of  the 
first  two  Important  coupled  vibration  modes  is  about  twice  as  great  for  the 
SWS-3-53  model  (average  g  of  0.04)  as  it  is  for  the  SWS-3b-53  model  (average  g 
of  0.02).  The  SWS-3  series  flutter  points  form  the  sharp  increases  in  the 
regions  of  instability  or  "buckets"  of  Figs.  3,  4  and  5;  thus,  the  mode  of  flutter 
may  be  one  that  Is  very  sensitive  to  g  variations.  Since  it  was  predicted  theo¬ 
retically  (Fig.  13)  that  the  mode  which  forms  the  "bucket”  Is  very  sensitive  to 
changes  in  g,  it  then  seems  possible  that  the  higher  structural  damping  of  the 
SWS-3-53  model  may  have  prevented  flutter  for  this  model  down  to  a  Mach 
number  of  1.8  where  It  was  destroyed  by  a  failure  of  the  inboard  leading  edge 
caused  by  a  root  seal  failure.  This  possibility  that  the  "buckets"  in  the 
experimental  curves  are  sensitive  to  g  variations  may  point  the  way  towards 
elimination  of  large  regions  of  instability  by  use  of  damping,  it  should  be  noted 
that  for  most  of  the  stabilizer  models  tested  the  value  of  g  for  the  first  two 
Important  coupled  modes  is  about  0.  02. 

The  data  for  the  SWS-3e-120  model  is  also  particularly  interesting  because 
this  model  failed  to  flutter  over  the  Mach  number  range  1.25  to  2,00.  This 

failure  to  flutter  means  that  the  curve  for  (w./w.  )2  =  0.  16  must  be  drawn  as 

N 

shown  in  Figs.  3.4  and  5  and  shows  that  at  these  higher  values  of  /w.  )  the 

0  "n 

bucket  is  not  evident.  Comparison  of  the  data  for  the  SWS-3e-120  and  the 
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other  SWS-3  models  helps  to  set  upper  and  lower  limits  of  (41^ 
in  the  Mach  number  range  1. 27  to  2. 10. 


/  <*>  ) 

aN 


for  flutter 


In  Fig.  6  and  Fig.  7  the  first  two  coupled  frequencies  with  the  pitching 
mechanism  in  operation,  and  w2,  were  used  to  form  the  flutter  parameters 
(b0.  75wi/af)  V<f*/*5)0  75  and  (bQ  75,  where  aj  is  the  speed 

of  sound  at  flutter.  The  use  of  the  relative  density  correction  in  this  form  is 
based  on  the  previous  experimental  results  of  Ref.  3  and  not  on  any  firm 
theoretical  basis.  Figures  6  and  7  may  be  useful  as  design  charts;  a  straight 
line  parallel  to  the  abscissa  being  a  constant  altitude  line,  and  a  straight  line  from 
the  suppressed  origin  being  a  line  of  constant  dynamic  pressure. 

In  Fig.  6,  the  first'coupled  vibration  frequency,  Wj,  is  used  to  normal  i  ze 
the  data.  This  vibration  mode,  as  can  be  seen  from  Table  5,  is  essentially 
a  combination  of  the  rigid  pitch  and  the  first  bending  modes  of  the  model.  The 
SWS-le-74  and  the  SWS-3d-87  both  have  the  same  value  of  the  parameter  (u>- /w^) 
and  hence  must  fall  along  the  same  boundary.  Thus,  the  curves  must  be  drawn 
as  shown  in  Fig.  8  with  a  narrow  stable  region  between  the  (o^/wj)  =  co  and  the 
1.60  curve. 

Figure  7  shows  curves  similar  to  those  of  Fig.  6  except  that  the  second 
coupled  vibration  frequency,  w2.  is  used  as  a  parameter.  This  vibration  mode, 
as  can  be  seen  from  the  data  of  Table  5  is  largely  a  combination  of  the  rigid 
pitch  and  first  torsion  modes  of  the  model  except  for  the  lowest  pitch  restraint 
stiffnesses  where  it  may  Involve  appreciable  bending. 

For  the  various  experimental  plots,  curves  have  been  drawn  on  the  basis  of 
a  bare  minimum  of  data.  The  fairing  of  such  curves  is  subject  to  some  question, 
and  Figs.  3, 4,  5, 6,  and  7  therefore  represent  only  rough  sketches  of  where  the 
flutter  boundaries  lie.  The  general  outlines  of  the  curves  are  probably  correct, 
and  enough  experimental  data  has  been  obtained  to  show  that  there  are  large 
Increases  in  the  regions  of  Instability  with  sufficiently  low  values  of  the  pitching 
frequency.  Furthermore,  these  increases  appear  to  follow  the  general  trends 
established  by  the  theoretical  results. 

In  one  respect  the  theoretical  results  do  not  match  the  experimental  results 
oven  qualitatively.  This  is  at  the  lower  Mach  number  of  the  calculation  M  =  1 .  52 
or  Mfl  =  10/9.  For  this  case  the  calculated  results  show  that  even  the  "locked" 
case  has  a  sharp  increase  in  the  region  of  instability  and  predicts  that  the  SWS-1 
and  the  SWS-3  series  models  will  flutter  in  the  cantilever  or  "locked"  configura¬ 
tion.  The  failure  of  the  theoretical  calculations  to  predict  flutter  correctly  in  this 
leglmc  Is  probably  due  to  the  failure  of  the  linearized  aerodynamic  theory  to 
predict  aerodynamic  forces  correctly  in  the  high-transonic  -  low  supersonic 
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regime.  Similarly,  it  seems  probable  that  the  failure  of  the  theoretical 
calculations  to  make  good  quantitative  predictions  throughout  the  Mach  number 
range  for  the  various  pitching  frequencies  and  model  stiffnesses  is  due  to  the 
failure  of  the  aerodynamic  terms  in  describing  accurately  the  actual  forces  on 
the  wing. 
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SECTION  III 
CONCLUSIONS 


Some  conclusions  may  be  drawn  from  the  theoretical  and  the  experimental 
results  of  the  present  program.  They  may  be  summarized  as  follows: 

1.  Three  basic  assumed  modes  appear  to  be  sufficient  to  define 
qualitatively  the  flutter  boundaries  when  the  velocity-component 
method  of  Ref.  8  Is  used,  These  modes  are  wing  first  bending, 
wing  first  torsion,  and  rigid  pitch.  Addition  of  wing  second  bend¬ 
ing  does  not  change  the  results  of  the  calculation  significantly. 

2.  For  low  margins  of  safety  in  bending-torsion  flutter,  the 
inclusion  of  a  high  frequency  pitch  mode  results  in  minor 
reductions  in  flutter  speed. 

3.  For  both  low  (7%  and  high  (45%  margins  of  safety  in  bending- 

torsion  flutter,  the  inclusion  of  a  critical  pitch  mode  =  1) 

causes  large  regions  of  instability  in  an  essentially  pitch¬ 
bending  flutter  mode  which  may  extend  as  high  as  M  =  2. 


4. 


5. 


The  theoretical  calculations  do  not  give  a  good  quantitative 
correlation  with  the  experimental  results.  The  theoretical 
calculations  predict  larger  regions  of  instability  than  are 
observed  experimentally.  They  also  predict  that  the  rapid  in¬ 
creases  In  the  regions  of  instability  will  occur  at  higher  values 
of  (wA/wu  )  than  were  observed  experimentally. 

v  n  j 


The  theoretical  calculations  indicate  that  the  mode  of  flutter 
which  causes  the  large  increases  In  the  region  of  instability 
may  be  very  sensitive  to  changes  in  structural  damping  coef¬ 
ficient.  Some  of  the  test  data  obtained  from  the  SWS-3  series 
models  confirm  this  conclusion. 
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APPENDIX  I 


THEORETICAL  CALCULATIONS 


l  Introduction 

In  setting  up  the  flutter  equations  for  the  all-movable  swept  stabilizer,  the 
authors  examined  the  relative  merits  of  the  strip-theory  method  (Ref.  7)  and 
the  velocity-component  method  (Ref.  8).  For  the  strip-theory  method,  the 
aerodynamic  forces  are  applied  to  sections  parallel  to  the  free-stream  while 
for  the  velocity-component  method,  they  are  applied  to  sections  normal  to  the 
elastic  axis.  The  former  method  is  more  rational  when  the  wing  ribs  are 
parallel  to  the  free  stream,  and  gives  a  better  representation  of  the  aero¬ 
dynamic  conditions  at  the  root  and  wing  tip.  The  latter  method,  however, 
appears  to  be  more  suitable  for  the  swept  stabilizer  model  which  derives  all  its 
stiffness  characteristics  from  a  single  spar.  The  simple  spar  type  of  construc¬ 
tion,  the  relatively  high  length  to  chord  ratio  as  well  as  the  results  of  vibra¬ 
tion  tests  suggest  that  the  concept  of  the  root  being  effectively  clamped  perpen¬ 
dicular  to  the  elastic  axis,  which  is  a  basic  assumption  of  the  velocity-component 
method,  is  well  justified.  Therefore,  it  was  decided  that  the  velocity-component 
method  would  be  used  in  deriving  the  equations  of  motion. 

In  the  derivation  and  solution  of  the  equations  of  motion  by  the  velocity  com¬ 
ponent  method  all  quantities,  mass  parameters  and  aerodynamic  forces,  are  re¬ 
ferred  to  a  reference  system  (x  Q ,  yfl  )  swept  with  the  elastic  axis  (Fig.  8). 

In  particular  the  Mach  number  used  in  obtaining  the  aerodynamic  coefficients 
must  be  the  crossflow  Mach  number  Mg  .  In  the  presentation  of  the  results, 
however,  all  the  theoretical  flutter  parameters  have  been  referred  to  an 
unswept  reference  system  (x,  y)  for  convenience  when  comparing  with  experi¬ 
mental  results. 

2.  Flutter  Equations  Based  on  Velocity-Component  Method 

The  flutter  equations  are  derived  following  the  method  of  Section  16.2  of 
Ref.  9,  The  assumption  that  the  wing  displacement  is  a  superposition  of  four 
modes  gives  as  the  deflection  of  any  point  (Fig.  8) 
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la (xa-  =  Fhj  +  Fh2  ^nyh2w  *  xoF«(  W° 


+  (yQ  sin  g  +Xq  cos q  -d)  0  (t) 
where  (see  Eqs.  40-42) 

F^  ,  are  the  first  and  second  assumed  cantilever  bending  modes 
1  2 


(1) 


F^  is  the  first  assumed  uncoupled  torsion  mode 

Hj,  Ej  )  are  reference  tip  amplitudes  for  the  first  and  second  bending  modes, 
a  q  ,0  I  first  uncoupled  torsion  mode,  and  rigid  body  pitch  mode,  respectively. 


/ 

''"V 

/ 


F i g.  8,  Axis  system  for  swept  stabiliser 

From  Eq.  (1)  it  is  seen  that  the  rigid  body  pitch  is  equivalent  to  a  bending  of 
the  elastic  axis  plus  a  rotation  about  the  elastic  axis,  so  that  only  the  aero¬ 
dynamic  forces  due  to  the  translation  and  rotation  of  sections  normal  to  the 
elastic  axis  are  needed.  Application  of  the  Lagrange  equations  of  motion  to  the 
system  as  given  by  Fig.  8  along  with  the  assumption  of  simple  harmonic  motion 
and  the  introduction  of  the  dimensionless  variable 


’to  B 


leads  to  the  following  dimensionless  set  of  flutter  equations: 


DO 


0 


(2) 


(3) 
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4  b  dF. 
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(23) 


hui  x  h 

(26) 

hlh- 

=  -ltanfl  L» 
k  ^ 

(27) 

h.« 

•  L„  -  H  <5  *  »> 

(28) 

hia' 

(29) 

M«h 

.  Mh-Lhli  .  .1 

(30) 

M«h' 

’  *  a)] 
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(32) 


-  (Mh  ♦  L  J(I  +  a)  ♦  Lh(I  ♦  a)2 
or  or  a  n  or  2  •'2 


M  .  =  -i  I  -  -i  — 

0,(1  k  8  2k 


-  -i- - Lh(i  -  a2)]* 


H*  * 


\q  sin  Q- d 


hih  4  \a  cosfl  4  hih*  sinfl 


(33) 

(34) 


M 


or* 


ya  sin  Q  -  d 


M 


arh  4  M«a  cos  O  ♦  Mwh,  sin  O 


(35) 


where  L^,  Lg,  Mfl  and  Mh  are  as  defined  by  Ref.  10. 

The  above  equations  were  written  using  the  actual  mass  distribution  for  the 
stabilizer  being  studied.  These  relationships  are 


where 


(37) 


(38) 


(39) 


However,  to  simplify  the  aerodynamic  calculations,  the  tapered  planform  was  re¬ 
placed  by  a  rectangular  planform  of  constant  chord  so  that  the  aerodynamic  coef¬ 
ficients  would  remain  constant  along  the  span  at  a  given  value  of  reduced  frequency. 
A  check  calculation  has  shown  that  if  reference  semichord,  bf,  is  taken  at  the  75% 
span  station  of  ihe  actual  mode  perpendicular  to  the  elastic  axis,  the  difference 
between  the  values  of  the  aerodynamic  integrals  as  given  by  the  rectangular 


•It  should  be  noted  that  this  equation  is  given  incorrectly  in  Ref.  9. 
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planform  and  those  values  found  by  an  actual  numerical  integration  along  the 
span  of  the  tapered  wing  are  very  small. 


The  first  bending  mode  was  taken  as 


(40) 


the  first  torsion  as 

Fo  =  ’to 


(41) 


and  second  bending  as 

Fh  =  -12. 209  rig  +  25.488rj*  -12.279tj^  (42) 

2 


The  second  bending  mode  was  obtained  by  assuming  a  power  series  in  t]q  which 
satisfied 

(1)  the  boundary  conditions  for  a  cantilever  mount, 

(2)  the  condition  of  orthogonality  with  the  first  bending  mode,  and 

(3)  the  condition  of  zero  deflection  at  the  75  percent  span  location. 

Condition  (3)  was  obtained  from  observation  of  the  node  line  for  the  second 
bending  mode  of  the  actual  mode  during  vibration  tests. 


The  parameters  used  in  the  analyses  were 

=  30  a  r  _0. 20 


r*  -  0.  250 
nQ 


S3 -  43°14' 


=  0. 20 


=  0.  66798 


0.21233 


sin  M 


3. 22603 
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which  resulted  in  the  following  values  for  the  coefficients  of  the  flutter  equations. 


A  =  0.078,125  ^  +  0.025,919  Lhh,  +  2.071,439  l-(-—  )  Z|  («) 

8  =  0.020,3321^  +  0.051,9581^,  (44) 

C  *  0.061,035  +  0.010,799  1^,  +  0,342,857  (45) 

D  =  -0.228,066  +  0.059,291  Lho  +  0.055,743  L^,  +  6.582,193  (46) 

E  =  0.020,331,530  Lhh  -  0.000,120,555  Lhh,  (47) 

F  =  0.091,120  +  0.025,919  Lhh,  +  3.417,742  1  -  J  2  Z  (48) 

G  = -0.000,283,89  Lhfl  -  0.004,972,8  1^,-0.212,126  (49) 

H  =  0.038,584  Lj^  -0.027,302  Lhfl  -0.025,668  Lhh,  -1.161,620  (50) 

I  =  0.061,035  Mflh  +  0.021,599  Mah,  +  0.  342,  857  (51) 

J  =  -0.000,283,89  Mflh  +  0.037,  371  Mah,  -0.  212, 126  (52) 

K  =  0.060, 863  M  +  0.010, 124  M  .  +  0.  441 , 964  (1  -  Z  j  (53) 

L  = -0. 180,995  Mflh  +  0.055,585  Moa  +  0. 052,  259  Mffh.  +  1.  691, 275  (54) 


M  =  0.059,291  Mah  +  0.023,604  Mah,  +  0.228,066  Lhh  +  0.076,860  Lhh>  +6.582,193 

(55) 

N  =  -0.027,302  Mflh  + 0.023,604  Mah.  +0.038,584  Lhh  +  0.158,271  Lhh,  -  1.161,620 

(56) 

O  =  0.055, 585M  +  0.014. 753M  ,  ♦  0.180,995  L._  +  0.  030,618  L.  .  ,  +  1.691 ,275 

aa  aa  nrr  nn 

(57) 

P  =  -0. 168,098  Moh  +0.080,996  Mffa  +  0.076,149  Mffh,  -  0.687,648  Lhf) 


♦  0. 168,098  Lha  *  0.  158,038  L^,  +  23.195.949  1  -^A  j2Z  j  (58) 


3  Solution  of  Equations  of  Motion  for  Flutter 

The  flutter  equations  for  the  swept  stabilizer  were  solved  for  two  (bending- 
torsion),  three  (bending-torsion -pitch)  and  four  (first  bending-second  bending- 
torsion-pitch)  degree -of -freedom  systems.  The  flutter  determinants  for  each 
system  are  respectively 
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where  A,  B, - P  are  given  by  Eqs.  (43)  through  (58).  In  each  case  the 

aerodynamic  terms  were  evaluated  by  selecting  specific  combinations  of  M  and 
k. 


The  general  pattern  of  solution  of  the  three  flutter  determinants  was  the 
same.  A  given  determinant  was  first  expanded  into  a  complex  polynomial.  Since 
the  right-hand  side  of  the  equation  was  zero,  two  separate  equations  were 
written  by  setting  both  the  real  and  Imaginary  parts  of  the  polynomial  equal  to 
zero.  These  two  simultaneous  equations  were  solved  for  any  two  desired 
eigenvalues. 


In  the  iwo-degree-of-freedom  case,  the  complex  polynomial  resulting  from 

2 

the  expansion  of  the  determinant  was  solved  for  the  eigenvalues  (wh  /w  )  and 

2  1  * 

(w  /u>j)  .  The  cross  flow  Mach  numbers  used  were  Mg  =  0,  10/9,  5/4  and 

10/7.  This  case  corresponds  to  the  pitch-locked  condition. 


The  three-degree-of-freedom  system  was  solved  for  the  two  eigenvalues 
2  2  2 

and  (i i) n  /u>j)  .  (w^  /wffl  )  was  set  equal  to  0.0625,  a  value  which 

corresponded  closely  to  the  average  value  for  the  actual  stabilizer  models. 

Again  Mg  ^  0,  10/9,  5/4,  and  10/7  were  used. 

2 

Finally,  for  the  four-degree-of-f reedom  systems,  Z  =  (w  /wf)  (1  +  ig)  was 

ON  I 

*  2 

used  as  the  eigenvalue.  Here,  it  was  necessary  to  specify  values  of  (uj.  /w  )  , 

n ,  O'  ■ 

2  2  11 
U h  /vt>  )  and  (u>A/a>  )  in  advance.  The  solutions  of  the  fourth-order  deter- 

"2  i  v  "i 

minants  were  carried  out  on  a  650  IBM  computer  using  a  program  developed  by 
North  American  Aviation  in  Columbus,  Ohio,  and  the  results  plotted  on  a 
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V-g  diagram.  .  Value  of  the  constants  for  which  solutions  were  found  were 

(«.  /V  )2  =  0.0625,  (on  /u>  )2  =  0.9025,  1.155625  for  (w./w  )2  =0.30 

nl  “1  n2  “l  v  al 

and  («.  /uj  )2  »  0.0625,  (on  /w  )2  =  1. 155625  for  (w./w  )2  =  0.20. 

l>i  OTj  "2  Q1  ♦  °l 

The  results  from  the  two-  and  three-degree-of-freedom  cases  are  plotted  in 

Fig.  9  and  then  crossplotted  in  Fig.  1.  The  pitch-locked  values  shown  as  as- 
2  2 

ymptotes  as  )  — »•  ae  In  Fig.  9  and  as  the  (w^/  )4— »  ®  boundary  in 

Fig.  1  are  the  results  of  the  two-degree-of-freedom  calculations. 

The  most  interesting  feature  of  these  analyses  are  the  very  deep  "buckets" 
that  occur  at  low  values  of  (w./w^  in  Fig.  9  and  correspondingly  at  low 
2  1 

values  of  (u>^/ ln  Fig.  1.  In  some  cases  the  curves  actually  double  back 

on  themselves  giving  two  regions  of  stability  at  a  given  value  of  (<*>,/ u>.  )2  or 

*  1 

Mach  number.  The  presence  of  these  "buckets"  is  apparently  due  to  a  change  in 
flutter  mode  shape  and  can  be  explained  by  looking  at  sample  four-degree-of- 
freedom  calculations  in  some  detail. 


Each  solution  for  the  four-degree-of-freedom  problem  at  a  given  set  of 

values  for  (w.  /«  )2,  (wh  /u>  )2,  and  (wA/w  )2  and  Mach  number,  yields 

nl  al  h2  °1  ♦  “l 

four  separate  curves  on  branches  on  a  V-g  plot  and  several  values  of  k  for  the 
flutter  condition  of  g  =  0  (see  Fig.  10).  Since  each  branch  represents  a 
particular  mode  of  flutter,  it  appears  from  Fig.  10  at  =  10/9  (M  =  1.  525) 

that  the  stabilizer  is  capable  of  flutter  in  the  1st  mode  and  the  2nd  mode.  At 
M  q  -  5/4  (M  a  1.  716)  the  stabilizer  has  one  unstable  region  along  the  V  axis 
in  the  2nd  mode  and  two  unstable  regions  in  the  3rd  mode.  A  set  of  three- 
dimensional  sketches  of  V  versus  g  versus  M  is  shown  in  Fig.  11.  To  avoid 
confusion,  each  sketch  contains  only  one  type  of  flutter  mode.  Because  of  the 
difficulty  in  following  the  various  possible  flutter  modes  from  a  V-g  diagram  to 
a  Vj/tu^  bg.  versus  M  plot,  the  results  of  the  four-degree-of-freedom 

analysis  of  Fig.  12  were  ultimately  drawn  after  looking  at  three-dimensional 
plots  of  V  versus  g  versus  M  with  interest  concentrated  on  the  traces  of  the 
different  modes  in  the  g  =  0  plane. 


The  lowest  set  of  curves  on  Vj/w^  bQ  ?5  versus  M  in  Fig.  12  form  the 

critical  flutter  boundary.  This  boundary,  as  can  be  seen  by  looking  at  Fig.  12 
is  formed  by  three  different  flutter  modes  each  becoming  the  critical  boundary 
of  instability  over  a  particular  Mach  number  range.  The  flutter  boundary  from 
the  four-degree-of-freedom  analysis  for  (u>0/wa  )2  =  0.  30  appears  to  compare 
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MODE  (2) 


© 

a 


© 

fH 

O 


Fig.  12.  Flufler  potometers  bQ  ^  ond  V^/V^Bq  versus  Mach  number  from  four-degree-of-freedom 

co^ulotions. 
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very  favorably  with  the  three-degree-of-freedom  flutter  boundary 


(tu.  /w  )‘ 


which  is  also  shown  in  Fig.  12.  Thus,  the  second-bending  degree  of  freedom 

Vf 

apparently  has  little  Influence  on  the  flutter  parameter  -  of  a 

<‘,n1b0.75 

swept  stabilizer  below  M  *  2. 0  except  for  a  general  lowering  of  the  curve.  In 
using  the  four-degree-of-freedom  analysis  to  interpret  the  points  found  in  the 
three-degree-of-freedom  analysis,  it  is  seen  that  the  8-shaped  curves  of  Fig. 
do  indeed  appear  reasonable  and  are  a  direct  result  of  a  change  in  critical 
flutter  modes  in  going  from  low  Mach  number  to  high  Mach  number. 


9 


Another  Interesting  characteristic  of  the  V-g  solutions  of  the  four-degree- 
of-freedom  analysis  is  the  variation  in  the  flutter  boundary  with  small  changes 
in  the  structural  damping  coefficient,  "g.  "  By  referring  to  Fig.  13,  it  is  seen 
that  increasing  the  structural  damping  from  g  *  0  to  g  =  0.  06  moves  the  "bucket” 
on  the  flutter  boundary  due  to  the  2nd  mode  from  about  M  »  1.75  to  about 
M  *  1. 65.  The  general  level  of  the  flutter  boundary  as  determined  by  the  1st 
and  2nd  modes  will  not  be  changed. 


Fig  13.  Vocation  of  Vp-'»,  bp  versus  Mach  number  with  changein  structural 
damping  from  four-degree-of  ^•wdom  calculations. 
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APPENDIX  II 


EXPERIMENTAL  DATA 


This  appendix  gives  the  detailed  tabulation  of  both  design  and  experimental 
data.  Since  both  the  model  design  and  the  testing  techniques  are  essentially  the 
same  as  those  described  in  Ref.  2,  little  discussion  of  them  is  included  in  this 

appendix. 

The  planform  of  the  stabilizer  models  is  shown  in  Fig.  14a  and  a  cross 
section  of  the  root  of  the  SWS-1  series  models  is  shown  in  Fig.  14b.  Root 
cross  sections  for  the  SWS-2  and  the  SWS-3  series  models  are  nut  shown  since 
they  differ  only  in  minor  details  from  the  SWS-1  model.  As  can  be  seen  in 
Fig.  14b  the  spar,  which  contributes  essentially  all  of  the  model  bending 
stiffness  and  most  of  the  torsional  stiffness,  is  constructed  of  a  pine  core 
around  which  is  wrapped  an  aluminum  skin.  Steel  caps  are  then  cemented  to  the 
spai .  Both  the  steel  and  aluminum  are  tapered  linearly  along  the  span  giving, 
when  combined  with  the  taper  of  the  height  and  width  of  spar,  the  required  fourth 
power  distribution  to  the  bending  and  torsional  rigidities.  El^  and  GJfl  : 


(h°  \  =  (  1  -  V/2  £ )  (G2) 

) 

o 


Balsa  wood  cemented  to  the  spar  was  used  to  give  the  aerodynamic  shape  re¬ 
quired,  and  suitably  spaced  lead  weights  were  used  to  give  the  mass  parameters 
required.  Table  1  gives  a  summary  of  the  design  parameters  for  all  cf  the 
stabilizer  models. 
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Tcble  1.  Design  parumetefs  fo<  -.  vtpl  slobilner  models. 

The  pararartvrs  presented  in  this  tsble  *re  common  ro  nil  the  models  built  in  this  t*ojtraxu. 


Geometric  Parameters 

Panel  aspect  ratio,  AR .  1-2  1 

Taper  ratio,  A . . .  1  2 

Sweep  angle  of  1/4  chord .  4V0° 

Mean  aerodynamic  chord  (in.),  UAC .  7.777P 

Section  mas.  thickness  (%  chord),  .  ('.O'; 

Line  of  max.  thickness  (%  chord) . . .  S0.CS 

Design  P ammeters 

Section  center  of  gravity  location  (%  chord),  (eg) .  50.0% 

Radius  of  gyration  (fraction  of  semichord),  r^. .  0.50 

Calculated  locus  of  shear  centers  (IS  chord),  (eo) .  40.0% 

Properties  of  Balsa  Wood  (average  values) 

Modulus  of  elasticity  in -bending  (lb/in2),  E .  6400 

Modulus  of  elasticity  in  torsion  (lb/in2),  G .  20,000 

Density  (lb/in^),  pBn .  0.00)900 

Properties  of  Pine  Core 

Modulus  of  elasticity  in  beuding  (lb/in2),  E .  1.429  -•  10 

Modulus  of  elasticity  in  tc.rion  (lb/in2),  G  .  0. 107  »  I0l 

Density  (lb/in3),  pp .  0.014 


The  root  fitting  and  the  mounting  block  with  the  pitching  mechanism  are 
shown  In  Fig.  15.  The  spar  was  glued  and  screwed  to  the  root  fitting,  shown 
removed  from  the  mounting  in  Fig.  15a.  Pitching  frequency  was  controlled 
by  changing  the  thickness  of  the  flexure  shown  on  the  end  of  the  root  fitting  in 
Fig.  15a.  Figure  15b  shows  the  rear  of  the  mounting  block  with  the  flexure 
In  place.  The  angle  of  attack  of  the  model  could  be  changed  by  rotating  the 
whole  clamp  shown  in  Fig.  15b.  Drag  and  lift  loads  were  carried  adequately 
by  three  ball  bearings  in  the  mounting  block.  The  gap  between  the  root  and  the 
mounting  block  was  sealed  with  aluminum  foil  for  all  tests. 

With  the  pitching  mechanism  "locked  Out,  "  static  tests  were  made  on  most 
of  the  models  in  an  attempt  to  determine  the  cantilever  properties  of  the  model. 
The  properties  determined  were  measured  elastic  axis,  (ea>M .  as  discussed  in 
Ref.  2.  and  the  root  values  of  El q  and  GJg.  The  results  of  these  measurements 
are  given  in  Table  2.  There  is  considerable  scatter  in  the  EI^  ,  GJ^  and  (ca)^ 
data. 


The  measured  mass  per  unit  length  at  the  root  (mo)M  is  also  given  in  Tabic  2, 

This  quantity  was  indirectlv  measured  using  the  assumed  mass  distribution 

2 

m  (y)  =  (mJM  (1  -  -i— .  )  (63) 

°  M  21 

By  just  measuring  the  total  mass  and  then  computing  (m  )..  from 

O  M 


total  mass 


f  1  (  1  -  -  )  2  dy 

o  21 


(64) 
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Fig.  15b. 


Fig.  IS.  Pi cluresol  root  mounting  block. 


40 


In  Eq.  (64)  the  total  mass  does  not  include  the  mass  of  the  root  fitting,  so  that 
the  value  for  (m0) j^j  includes  only  the  mass  of  the  balsa,  lead,  glue  and  the 

spar. 

Data  for  the  pitching  frequency  is  found  in  Table  3.  The  mass  moment  of 
inertia  of  the  whole  model,  including  the  root  fitting,  was  obtained  by  swinging 
the  model  with  a  bifilar  pendulum.  The  pitching  mechanism  flexibility  influence 
coefficients,  C^,  was  measured  with  a  transit  and  mirror  arrangement.  The 
pitching  frequency  was  then  calculated  as: 

f  =  -1-  J - L  ■  (65) 

The  results  of  the  flutter  tests  are  given  in  Table  4.  For  the  sake  of 
convenience,  most  of  the  important  experimental  natural  still-air-vibration 
frequencies  are  included  as  well  as  the  tunnel  conditions  at  flutter.  If  flutter 
occurred,  the  conditions  at  the  start  of  flutter  are  given.  If  no  flutter  occurred 
during  the  test  run,  the  conditions  at  the  start  and  end  of  the  test  are  given. 
Figures  16  and  17  are  excerpts  from  the  high  speed  movies  taken  during  the 
flutter  of  the  SWS-1-98  and  the  SWS-3d-8?  models,  respectively.  These  portions 
of  the  movies  have  been  analyzed  and  the  results  are  presented  in  terms  of  the 
pitching  motion  at  the  root  and  the  motion  of  the  Up  sections  in  Figs.  18  and  19. 
These  flutter  modes  are  typical  of  those  encountered  for  the  stabilizer  models. 

Complete  vibration  data,  including  sketches  of  node  lines,  frequency,  and 
structural  damping  of  the  lower  modes  of  vibration  are  found  in  Table  5.  All 
of  the  models  were  vibration  tested  in  both  the  "locked,  "  or  cantilever,  condi¬ 
tion  and  with  the  pitching  mechanism  in.  Figure  20  is  a  plot  of  the  normalized 
coupled  frequencies  with  the  pitching  mechanism  in.  The  lowest  cantilever 
bending  frequency,  fj^,  was  used  as  the  normalizing  frequency.  It  is  Interest!^ 
to  note  that  the  frequencies  of  the  first  coupled  modes,  fj,  for  most  of  the 
stabilizer  models  fall  along  a  common  curve  with  not  too  much  scatter.  The 
same  is  true  for  the  frequencies  of  the  second  coupled  mode,  fg. 

Table  6  gives  the  influence  coefficient  data  for  the  models  with  the  pitching 
mechanism  in  and  Fig.  21  shows  the  location  of  the  stations  at  which  influence 
coefficients  were  taken. 
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Tobl*  2.  Static  data  for  swap)  stabilizer  models. 


Model 

K>m 

(*luj/ft) 

(ee}M 
(*  chord) 

<CJon)M 

(E'^ 

STS-1 

0.0214 

3.690  y  104 

8.493  y  104 

STS- 2 

0.0214 

2.696  y  104* 

6.860  y  104* 

STS- lb 

0.0220 

2.877  >:  104 

6.661  y  104 

STS-lc 

0.0223 

46.3* 

3.459  y  104 

7.422  y  104 

STS- Id 

0.0225 

49.0* 

2.380  y  104 

5.84  y  104  ! 

STS-  le 

0.0246 

39.0* 

4.220  *  104 

7.836  y  104 

STS-} 

0.0240 

33.36* 

7.572  y  104 

10.448  y  104 

STf  }• 

0.0231 

43.0* 

3.680  y  104 

8.59  y  104 

STS-  3b 

0.0224 

50* 

5.23  y  104 

7.139  y  104 

STS- 3c 

0.0233 

1.942  y  104* 

6.754  y  104* 

STS- 3d 

0.0247 

50* 

5.234  y  104 

7.489  y  104 

STS-Je 

0.0274 

3-89  y  104* 

8.89  x  104» 

*  Dui 

for  spar  only. 

Table  3.  Pitching  frequency  data. 


Model 

(•lu**ft2) 

(red/lb-h) 

»* 

(cp») 

STS-l-L 

0.00244 

0.000022 

5TS- 1-138 

0.00244 

0.000343 

138 

STS- 1-105 

0.00244 

0.000915 

103 

STS- 1-98 

0.00244 

0.001072 

98.3 

STS-2-L 

0.000022 

STS- lb- L' 

0.00223 

0.000022 

STS-lc-48 

0,00238 

0.004460 

47.6 

STS-ld-L 

0.0000>2 

STS-le-74 

0.00226 

0.002132 

'3.7 

STS- 3-5  3 

0.00206 

0.004460 

52.5 

STS-3.-63 

0.00224 

0.002854 

63.1 

STS-}b-53 

0.00212 

0.004249 

53.0 

STS-}c-74 

0.00205 

0.00224 

74.4 

ST5-W-8? 

0.00225 

0.001498 

86.7 

STS*  5e- 120 

0.00 275 

0.000628 

120.0 

WADC  TR  56-285 


42 


Table  4.  Experimental  flutter  data. 
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Fig  ?0  Vibration  frequency  dato  for  swept  stobilirer  rrodels 
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Fij.  ?1.  Localion  of  influence  coefficient  notions. 
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Table  5.  Expcrimanfal  vibration  data. 


Toble  5.  Experimental  .vibration  dofq.  (cont.) 
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Tobl*  6.  E*p*f influence  coefficient  data.  SWS- !W- 142 
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